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THE  MORPHOLOGY  OF  A  MULTI-BUBBLE  SYSTEM 
IN  THE  IONOSPHERE 


I.  INTRODUCTION 

The  behavior  of  Ionospheric  plasma  has  been  extensively  studied  with 
the  aim  of  understanding  various  ionospheric  phenomena  including  equatorial 
spread  F  (ESF)  and  plasma  cloud  strlations  (see  the  following  reviews: 
Ossakow,  1979;  Fejer  and  Kelley,  1980;  Ossakow,  1981;  Kelley  and  McClure, 
1981;  Ossakow  et  al.,  1982).  In  particular,  ESF  is  thought  to  be  initiated 
by  the  Rayleigh-Taylor  instability,  first  proposed  by  Dungey  (1956).  In 
the  context  of  this  idea,  the  plasma  density  depletions  resulting  from  the 
instability  acquire  upward  polarization  induced  E_  x  B_ drift  velocities.  A 
large  body  of  literature  has  since  been  developed  to  describe  the  linear 
and  nonlinear  properties  of  Rayleigh-Taylor  plasma  density  irregularities 
under  various  assumptions  [Haerendel,  1975;  paisley  et  al.,  1972; 
Chaturvedi  and  Kaw,  1975a, b;  Hudson  and  Kennel,  1975.].  In  particular, 
considerable  attention  has  been  given  to  the  morphology  and  motion 
[Scannapleco  and  Ossakow,  1976;  Ossakow  and  Chaturvedi,  1978;  Ott,  1978; 
Hudson,  1978;  Anderson  and  Haerendel,  1979;  Ossakow  et  al.,  1979;  Zalesak 
and  Ossakow,  1980;  Zalesak  et  al.,  1982]  of  plasma  depletions  ("bubbles") 
in  the  equatorial  ionosphere.  Moreover,  a  number  of  observations  [Woodman 
and  LaHoz,  1976;  Kelley  et  al.,  1976;  Hanson  and  Sanatani,  1971;  McClure  et 
al. ,  1977;  Szuszczewicz  et  al.,  1980,  1981]  have  indicated  the  presence  of 
rising  plasma  bubbles. 

Another  phenomenon  of  Interest  is  that  of  the  strlations  in 
(artificial)  plasma  "clouds"  (density  enhancements).  This  effect  has  been 
attributed  to  the  x  gradient  drift  instability  [Linson  and  Workman, 
1970;  Volk  and  Haerendel,  1971]  and  appears  to  be  amenable  to  treatments 
similar  to  ESF  [Scannapleco  and  Ossakow,  1976;  Scannapleco  et  al..  1976; 
Ossakow  and  Chaturvedi,  1978].  The  Rayleigh-Taylor  instability  and  the  E_  x 
J5_  gradient  drift  Instability  are  both  interchange  modes  that  may  occur  in 
the  leading  edge  of  plasma  bubbles  and  backside  of  clouds  respectively.  In 
both  cases,  the  resulting  density  depletions  and  enhancements  are  thought 
to  drift  by  the  polarization  induced  E_  x  B_  drift.  The  electric  field  is 
produced  by  the  polarization  of  the  plasma  across  the  earth's  magnetic 
field.  The  essential  ingredient  Ls  the  small  but  nonzero  ion-neutral 
collision  frequency. 

Manuscript  approved  November  24, 1982. 
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As  a  result  of  the  Initial  Instability  such  as  the  Rayleigh-Taylor 
instability,  an  array  of  density  depletions  are  formed  with  the  wave  "‘tctor 
perpendicular  to  the  earth's  magnetic  field.  Moreover,  the  leading  edge  of 
a  bubble  itself  is  Rayleigh-Taylor  unstable  resulting  in  further 
"bifurcation"  [Ossakow  and  Chaturvedl,  1978].  Similarly,  the  x  _B 
instability  causes  the  backside  of  an  initial  cloud  to  striate,  forming 
"finger-like"  structures.  Thus,  the  plasma  depletions  and  enhancements 
typically  occur  in  multitudes.  Indeed,  McClure  et  al.,  (1977)  and 
Szuazczewlcz  et  al«,  (1980; 1981)  have  given  evidence  for  multiple  bubbles. 
In  the  previous  theoretical  and  numerical  simulation  works  [Ossakow  and 
Chaturvedl,  1978;  Anderson  and  Haerendel,  1979;  Zalesak  and  Ossakow,  1980; 
Overman  et  al.,  1982],  the  morphology  and  evolution  have  been  studied  using 
one-bubble  models  or  a  uniformly  distributed  array  of  bubbles  (clouds).  As 
a  result,  the  intrinsic  influence  of  neighboring  bubbles  on  each  other  has 
not  been  well  quantified.  Because  the  bubbles  rise  due  to  the  polarization 
Induced  E_  x  B_  drift,  the  electric  field  configuration  in  and  around  the 
bubbles  is  of  central  importance. 

In  the  present  paper,  we  study  the  structure  of  electrically 
interacting  multi-bubble  systems  and  seek  to  Identify  the  nature  and 
effects  of  the  mutual  interaction.  For  this  purpose,  we  start  with  a 
simple  two-bubble  configuration  in  the  context  of  a  fluid  description  and 
with  the  emphasis  on  identifying  the  basic  physics.  A  simple  plecewis« 
constant  density  profile  is  used  and  the  electric  field  inside  and  outside 
the  bubbles  is  obtained  analytically  by  the  method  of  images.  The  method 
is  then  applied  to  multi-bubble  configurations.  In  a  sense,  the  present 
work  is  a  non-trivlal  extension  of  the  single  bubble  work  of  Ossakow  and 
Chaturvedl  (1978). 

The  scope  of  the  present  paper  is  limited  to  discussing  the  morphology 
of  multi-bubble  systems  and  the  time-dependent  evolution  is  not  explicitly 
considered.  However,  the  implications  of  the  results  will  be  discussed  in 
the  context  of  the  behavior  of  bubbles  consistent  with  the  approximations 
used  in  the  analysis.  Although  we  primarily  discuss  the  ESF  plasma  density 
depletions  (bubbles),  the  technique  developed  here  can  be  extended 
straightforwardly  to  the  treatment  of  plasma  density  enhancements  (clouds). 


In  Section  II,  we  develop  a  two-^bubble  model  based  on  an  electrostatic 
(dielectric)  analogy  and  solve  the  current  conservation  equation  with 
quasi-neutrality  for  the  electric  field.  Section  III  describes  a  multi- 
bubble  model  and  section  IV  contains  the  discussion. 

II.  A  TWO-BUBBLE  MODEL 

A.  Formulation 

In  the  present  paper,  we  consider  the  electric  field  configuration  of 
electrically  interacting  multi-bubble  systems  imbedded  In  a  uniform 
background  plasma  and  neutral  gas.  In  order  to  illustrate  the  basic 
physics  and  the  theoretical  method,  we  first  develop  a  simple  two-bubble 
model.  For  his  purpose,  we  adopt  a  sharp-boundary  density  profile  in  which 
the  plasma  density  is  piecewise  constant,  being  uniform  inside  (n^)  and 
outside  (no)  the  bubbles  and  having  a  discontinuity  at  the  bubble 
boundaries  [Haerendel,  1973;  Ossakow  and  Chaturvedt,  1978;  Overman  et  al., 
1982].  The  bubbles  are  modelled  by  two-dimensional  cylinders  at  the  same 
altitude  with  circular  cross-sections  and  the  axis  of  the  cylinders  arc 
aligned  with  the  earth's  magnetic  field  which  is  assumed  to  be  uniform 
along  the  positive  z-axis.  No  neutral  wind  is  included.  Figure  1  shows 
schematically  the  geometry  and  the  coordinate  system.  The  bubbles  are 
located  at  x  *  -  *c  and  x  ■  xQ  so  that  the  inter-bubble  separation  distance 
is  2xQ.  For  reference  purposes,  we  denote  the  cylinders  by  Cl  and  C2, 
respectively.  The  radius  of  each  bubble  is  a.  In  the  equatorial  and  low 
latitude  regions,  the  x-axis  is  along  the  east-west  direction  and  the 
gravitational  force  is  along  the  negative  y-axis. 

The  basic  equations  in  our  model  are  the  particle  conservation, 
momentum  conservation  and  current  conservation  equations  [see,  for  example, 
Ossakow,  1981].  In  addition,  quasi-neutrality  can  generally  be  assumed  on 
the  time  scales  of  interest.  In  the  present  paper,  we  study  the  morphology 
of  two-  and  multi-bubble  systems  by  solving  the  current  conservation 
equation  to  obtain  the  instantaneous  electric  field  perpendicular  to  the 
magnetic  field. 
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The  current  conservation  equation  with  quasi-neutrality  is 

7  •  J  -  0  (l) 

where  £  Is  the  plasma  current  density  due  to  Ion  and  electron  drifts. 
Neglecting  the  inertia  terms  In  the  momentum  equations  for  cold  Ion  and 
electron  fluids,  the  current  density  In  the  lab  frame  can  be  expressed  as 
(see,  for  example,  Ossakow  et  al.,  1979) 

a 

where  ■  eB/m^c  Is  the  Ion  cyclotron  frequency,  z_  Is  the  unit  vector 

along  the  earth's  magnetic  field,  £  Is  the  gravitational  acceleration,  n  Is 
the  plasma  density,  and  Is  the  lon-neutral  collision  frequency.  In 
(2),  we  have  used  the  fact  that  me/m^  <<  1  and  have  neglected  the 

electron  j£  x  z_  contribution.  In  addition.  In  arriving  at  (2),  we  have  made 
the  approximation  vin/^1  «  !•  In  the  F  region,  v^^/Q.  Is  typically  of  the 
order  of  10"^  or  less. 

The  second  term  In  the  square  brackets  In  (2)  gives  the  force- 

field  (£  and  E)  aligned  drift  currents  due  to  the  finite  lon-neutral 

collisions.  It  is  convenient  to  separate  the  electric  field  £  according  to 

~  mi 

£ 5  £  -  -r  * 

The  term  (-m^/e)  £  is  the  component  of  E  cancelling  the  drift  along  the 
gravitational  field  so  that  the  net  drift  perpendicular  to  the  magnetic 
field  Is  described  by  E_.  Then,  the  current  J  can  be  written  as 


±-  °5j 


where 


®  s  vm 


5  l  +  77  “  &  x  *• 


t 


’  }  I  .*•** 


Here,  ]}  and  4  *  -  g£  are  assumed  to  be  uniform.  Physically,  Ej  can  be 
thought  of  as  the  electric  field  driving  the  current  J_  perpendicular  to  the 
magnetic  field  in  the  frame  moving  with  the  velocity  »  -sil^  &.  *  L 
relative  to  the  lab  frame  and  0  may  be  identified  as  the  Pedersen 
conductivity  due  to  v  *  0. 

Equation  (1)  can  now  be  written  in  the  equivalent  form 


7  •  (nE)  -  -  —  (&  x  z)  •  7n.  (5) 

in 

Perkins  et  al.  (1973)  obtained  this  expression  and  noted  that  (5)  describes 
a  dielectric  Immersed  in  a  uniform  electric  field  E  where 


(6) 

in 

In  an  earlier  work,  Longmlre  (1970)  utilized  a  similar  magnetostatic 
anology  to  treat  the  motion  of  Isolated  ion  clouds.  Ossakow  and  Chaturvedl 
(1978)  used  (5)  to  analytically  study  a  single  bubble  system.  In  this 
dielectric  analogy,  _E  is  the  polarization  (self)  electric  field  of  the 
bubbles  in  the  uniform  field  Eq  and  E^  corresponds  to  the  total  electric 
field  (E^j  +  E)  satisfying  the  boundary  conditions  across  the  bubble 
boundaries. 


(oE  ),  ■  continuous 

(7) 

(Ej)j  •  continuous 


and  at  infinity  (x,y  ♦  ") 


Sj*V  (8) 

Note  that  we  have  implicitly  chosen  a  reference  frame  in  which  the  electric 
field  of  the  distant  undisturbed  ionosphere  is  The  symbols 

I  and  1  refer  to  the  directions  parallel  and  perpendicular  to  the  boundary 
surfaces,  respectively.  In  the  present  paper,  we  also  adopt  the  dielectric 
analogy  and  solve  the  current  conservtlon  equation  (1)  subject  to  the  above 
boundary  conditions  (7)  and  (8).  For  this  purpose,  it  is  illuminating  to 
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rewrite  equation  (1)  as 


V* (oE  )  -  0  (9) 

—  — J 

In  the  following  section,  we  describe  the  method  of  Image  dipoles  used  to 
solve  this  "dielectric  equation".  As  a  matter  of  notation,  Ln  the 
remainder  of  the  paper,  we  use  E_  without  the  subscript  J  to  denote  the 
solution  of  (9). 

B.  The  Method  of  Image  Dipoles 

The  problem  of  solving  Poisson's  equation  (9)  with  multiple 
disconnected  boundaries  Is  generally  difficult.  However,  In  the  case 
treated  here  with  circular  cross-sections,  the  dielectric  analogy  allows  us 
to  construct  an  exact  solution.  Consider  first  a  single  dlelectrLc 
cylinder  of  radius  a  centered  at  x  *  0  and  a  line  charge  density  q  located 
at  x  -  b  (|b|>a).  It  is  well  known  (Smythe,  1968)  that  the  Induced 
electric  field  outside  the  cylinder  Is  that  of  a  line  charge  q'« 

-  q(l-K)/(l+K)  located  at  *  *  a  /b  and  a  line  charge  -q'  located  at  x»0. 
The  quantity  K  is  the  ratio  of  the  dielectric  constant  (o^)  Inside  the 
cylinder  to  that  outsLde  the  cylinder  (o^ ) 


The  Induced  electric  field  Inside  the  cylinder  Is  that  of  a  single  line 
charge  q"  *  2q/(i-HC)  located  at  x  -  b.  If  we  replace  the  line  charge  q  by 

A 

a  line  dipole  moment  P^  *  P^x  which  Is  equivalent  to  two  equal  and  opposite 
line  charges  separated  by  a  vanishingly  small  distance,  then  we  find  that 
the  electric  field  due  to  polarization  of  the  cylinder  Is  that  of  a  single 
image  dipole  IP  given  by 

-  “  “  ^T+k)  71  £o*  (U) 

o 

located  at  x  -  az/b.  Note  that  no  Image  dipole  Is  present  on  the  axis  (to 

be  contrasted  with  the  line  charge  case)  and  that  P  and  P  are  collnear. 

— o  — 

pointing  In  the  opposite  directions.  The  Induced  electric  field  Inside  the 
cylinder  is  that  due  to  a  dipole  moment  £*  given  by 
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located  at  x  *  b. 


(12) 


We  now  consider  two  identical  dielectric  cylinders  of  radius  a 
centered  at  x  *  -  xQ  and  x  ■  xQ  (see  figure  1)  immersed  in  a  uniform 
electric  field  E^.  Suppose,  for  the  moment,  that  the  two  cylinders  are 
non-interacting •  Then,  for  the  purpose  of  calculating  the  polarization 
electric  field  outside  the  cylinders,  each  cylinder  may  be  replaced  by  a 
dipole  moment  located  at  x  -  -x0  and  x  •  xQ,  where 

1  1  _ir  o 

to  -  7  E„  i’  <13> 

and  K  in  defined  by  equation  (10).  The  components  of  the  self  electric 
field  are 


Ex  ■  2p0Lf(x+vy)  +  f(x*x0»y^» 


(i4) 


and 


Ey  "  2P0[h(x-Hco,y)  +  h(x-xQ,y)J, 


(15) 


where 


f(x,y) 


.  x2  -  y2 
(x2+y2)2 


and 


(16) 


h(x,y) 


2  x  y 
(x2+y2)2 


(17) 


Here  the  x-  and  y-components  of  the  electric  field  due  to  a  dipole  PQ 
at  x  -  xQ  is  2P0f(x-x0,y)  and  2PQh(x-x0,y),  respectively.  The  electric 
field  E*  Inside  the  bubbles  is 


E* 


1+K 


(18) 


In  the  remainder  of  the  paper,  asterisks  will  be  used  to  denote  the 
electric  fields  Inside  the  bubbles  and  dipole  moments  producing  the  fields. 
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We  now  allow  tfo*  two  cylinders  to  Interact  with  each  other.  In 
addition  to  the  uniform  external  field  E  ,  each  cylinder  experiences  the 
dipole  field  of  the  other.  As  a  result,  each  cylinder  Is  further  polarized 
giving  rise  to  linage  dipoles  as  given  by  (11)  and  (12).  Iterating  the 
method  of  linage  dipoles  described  above,  it  Is  straightforward  to  show  that 
the  total  field  outside  the  bubbles  Is  given  by 

E  -  -  E  +  1L  2P  [f(x+x  ,y )  +  f(x-x  ,y)l,  (19) 
x  o  n*0  nL  v  n,,7J 


and 


Ey  *n-0  2pJh(x'*'xn,y)  +  h<x_xn»y>l» 


(20) 


where  f  and  h  are  defined  by  (16)  and  (17).  Here,  for  n  *  0, 


_  .  r 1-K\  a2  _ 

Pn  =  '  W  Pn-1  * 
n 


b  =  x  +  x  , , 
n  o  n-1 * 


(21) 

(22) 


and 

a2 

x  -  - - — -  . 

n  o  x  +x 

o  n-1 

For  n  «  0,  we  have  xn  »  xQ  and  PQ  Is  given  by  (10). 

Similarly,  the  total  electric  field  inside  each  cylinder, 
located  at  x  •  xQ,  Is  found  to  be 


(23) 


say,  (C2) 


E* 

x 


+  I+K  nl0  2V(X+Vy)' 


(24) 


and 


EJ  ■  -Rk  „Io  2p„h<x^n>y)  .  (25) 

where  Pn  and  ^  are  defined  above.  For  the  other  bubble  (Cl),  the  field  Is 
obtained  by  replacing  x0  with  -xQ  In  the  functions  f  and  h. 
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Ic  can  be  seen  from  (19)-(23)  that  the  convergence  properties  of  these 
series  depend  on  the  parameter  s  defined  by 


-  ilK 

8  =  1+K 


(2xoy 


(26) 


Each  successive  line  dipole  moment  Is  reduced  by  a  factor  of  s  and  a 
geometrical  factor  of  order  unity  from  the  preceding  one.  Since  xQ  >  a,  we 
have  Is |  <  1/4  for  any  K.  As  a  result,  the  series  are  generally  rapidly 
convergent  except  for  very  small  center-to-center  separation  distances 
(xQ  ~  a). 

C.  Applications 


In  order  to  apply  the  above  dielectric  results  to  the  Ionospheric 
bubble  and  cloud  problems  In  accordance  with  the  dielectric  analogy 
(equations  (5)  and  (6)),  we  identify  the  dielectric  constant  o  with  the 
Pedersen  conductivity  defined  In  Section  II  and  recall  that  E^  is  given  by 
equation  (6).  Note  that  the  Pedersen  conductivity  is  proportional  to  the 
plasma  density  so  that 


where  n^  and  ^  are  the  plasma  densities  inside  and  outside  the  bubbles. 

With  these  Identifications,  equations  (19),  (20),  (24)  and  (25)  describe 

the  electric  field  of  a  two-bubble  system  In  the  frame  moving  with  V,  ■ 
-1 

"  a  i  4  *  £  relative  to  the  earth.  Note  that  our  formalism  guarantees  that 
the  field  components  obtained  above  exactly  satisfy  equation  (9)  and  the 
boundary  conditions  (7)  and  (8),  as  can  be  verified  easily.  Thus,  the 
solution  Is  unique.  It  is  also  worth  noting  that  all  the  image  line  dipole 
moments  are  aligned  with  the  x-axls  and  no  higher  multipoles  such  as 
quadrupole  moment  arise.  If  the  bubbles  are  at  different  heights,  the 
Interaction  can  still  be  expressed  as  a  series  of  Image  dipole  moments  that 
are  parallel  to  the  x-axls  but  the  Images  are  not  Induced  on  the  x-axis 
Itself. 

As  a  result  of  the  polarization  electric  field,  the  plasma  bubbles 
(density  depletions)  E  x  B>  drift.  The  drift  velocity  relative  to  the 
distant  undisturbed  Ionosphere  where  the  electric  field  Is  Eg  is  given 
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2 

by  V  -  c  (E*  -  E  )  x  B/B  .  In  particular,  the  single-bubble  rise 
—  —  — o  — 

velocity  is  (using  equation  18) 


For  bubbles,  K  <  1  and  is  upward.  This  result  has  been  obtained  by 
Osgakow  and  Chaturvedl  [1978].  Because  £  is  assumed  to  be  uniform,  a 
single  isolated  bubble  maintains  its  circular  cross-section  as  it  rises 
with  the  constant  velocity  .  The  presence  of  a  second  bubble,  however, 
modifies  the  rise  velocity  significantly.  In  particular,  it  is  no  longer 
uniform  and  the  cross-sections  do  not  remain  unchanged.  We  have 
numerically  carried  out  the  summations  indicated  in  (19),  (20),  (24),  and 
(25).  The  results  are  depicted  in  figure  2.  The  solid  lines  represent  the 
self  electric  field  lines  or  equivalently  the  Pedersen  current  lines  while 

A 

the  dashed  lines  correspond  to  the  Instantaneous  E  x  J1  drift 
velocity  (Vg)  for  a  two-bubble  system  (2xQ  ■  2.5a  and  K  ■  0  for  100% 
depletion).  This  figure  shows  only  one  quadrant;  the  actual  system  is 
spatially  uniform  in  the  z-dlrection  and  is  symmetric  about  the  y-z  plane 
and  the  x-z  plane.  The  solid  electric  field  lines  are  such  that  the  line 
density  is  proportional  to  the  field  strength.  This  figure  clearly  shows 
that  the  electric  field  inside  the  bubble  is  significantly  modified  from 
the  uniform  field  of  an  isolated  bubble  (see  (18)).  As  a  general  remark, 
the  interaction  vanishes  as  K  >  1  and  the  separation  distance  increases  to 
infinity. 

In  this  example,  the  electric  field  strength  at  the  point  A  (x  *  x0~a, 
y  ■  0)  is  approximately  1/3  of  that  at  B  (x  »  xQ+a,  y  *  0)*  The 
prominently  nonuniform  electric  field  inside  the  bubbles  has  a  number  of 
Important  Implications  for  the  dynamic  behavior  of  the  two-bubble  system. 
The  drift  lines  (dashed)  in  figure  2  show  that  different  regions  of  a 

A 

bubble  undergo  drift  in  the  east-west  (x)  direction  with  respect  to  the 
undisturbed  ionosphere.  This  horizontal  drift  can  be  a  significant 
fraction  of  the  vertical  rise  velocity  at  some  points  inside  the  bubble. 
For  example,  |E*/(E*-Eq)|  -  0.55  at  r  -  a_,  and  6  -  ir/8  where  r  is  the 
radial  distance  from  x«x0,  a_  is  just  Inside  the  bubble  surface  and  8  is 
measured  from  the  point  A  (0  -  0).  Here,  xQ  -  1.25  and  K  »  0  have  been 
used  (figure  2).  The  reason  for  the  strong  divergence  in  the  field  lines 
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A  drawing  of  the  electric  field  lines  (solid  lines)  and  the  x  B_ 
drift  lines  (dashed  lines)  showing  one  quadrant  of  a  two-bubble 
system.  Tne  system  is  symmetric  about  the  x-z  plane  and  y-z 
plane.  The  separation  distance  Is  2xQ  ■  2.5a  and  K  -  0  (100X 
depletion).  The  points  A  and  B  are  Inside  the  bubbles  at  *  ■  xQ-a 
and  x  *  xQ+a,  respectively. 


Is  that  the  dipole  field  due  to  one  bubble  opposes  the  internal  field  of 
the  other  bubble* 


Figure  3  shows  the  two-bubble  rise  velocity  (Vy)  at  a  number  of  points 
inside  the  bubble  (C2)  relative  to  the  uniform  rise  velocity  Vj^  of  a  single 
bubble  given  by  (27).  The  ratio  R2  =  Vy/V^  is  plotted  for  the  interior 
points  A  and  B  (figure  2)  as  a  function  of  xQ/a.  We  note  that  the  rise 
velocity  Vy  is  strongly  affected  by  the  neighboring  bubble  for  center-to- 
center  separation  distances  (2xQ)  smaller  than  5a  to  6a.  Moreover, 
comparing  Vy  at  A  and  at  B,  we  see  that  the  E  x  JJ  drift  velocity  is 
sheared.  This,  together  with  the  conclusions  of  the  preceding  paragraph, 
implies  that  an  initially  circular  cross-section  would  not  remain  circular 
so  that  the  two-bubble  system  described  here  does  not  correspond  to  a 
steady-state  system.  This  should  be  contrasted  with  single-bubble  models 
in  which  steady-state  solutions  are  possible  for  piecewise  constant  density 
profiles.  Because  a  neighboring  bubble  generally  introduces  nonuniformity 
in  the  electric  field,  it  seems  difficult  to  construct  a  steady-state  two- 
bubble  system  unless  xQ/a  is  large. 

In  interpreting  the  results,  we  note  that  our  results  do  not  include 
the  time-dependence.  The  above  field  configuration  exists  if  the  two- 
bubble  system  as  described  is  created  at  some  time  (say,  t  ■  0).  Thus,  it 
would  be  appropriate  as  a  consistent  initial  condition  for  the  purpose  of 
studying  the  subsequent  time-evolution.  In  practice,  the  drift  lines  shown 
in  figure  2  are  expected  to  closely  approximate  the  actual  evolution  for 
some  period  of  time  after  t  ■  0  until  the  distortion  changes  the  field 
topology  significantly.  However,  the  general  feature  of  the  field  with  the 
weakest  x-component  and  hence  the  slowest  rise  velocity  in  the  region 
nearest  to  the  neighboring  bubble  should  remain  unchanged  in  time. 


III.  A  MULTI-BUBBLE  MODEL 

In  the  preceding  section,  we  have  discussed  in  detail  a  two-bubble 
model.  In  considering  a  multi-bubble  system  that  may  be  applicable  to  the 
ESF  phenomenon,  the  basic  physics  and  the  theoretical  treatment  remain 
unchanged.  However,  as  the  number  of  bubbles  Increases,  so  does  the  number 
of  image  dipoles.  For  an  N-bubble  system,  the  nc^  order  expression  must 
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Fig.  3  The  ratio  R2  of  the  vertical  rise  velocity  Vy  of  a  two-bubble 
system  to  the  single-bubble  rise  velocity  (equation  (27))  as  a 

function  of  xQ/a.  The  separation  distance  Is  2xQ  “  2.5a  and  K  “  0 
(100%  depletion).  The  curve  A  corresponds  to  the  point  A  (x  ■  xQ- 
a)  and  the  curve  B  to  the  point  B(x  ■  Xq+s)  In  figure  2. 
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Include  N(N-l)n  image  dipoles.  Fortunately,  the  Influence  of  a  bubble 

decreases  as  the  inverse  square  of  the  separation  distance.  Thus,  the 

nearest  and  the  second  nearest  neighbors  are  expected  to  have  the  dominant 

effects.  A  consideration  of  equations  (11)  and  (12)  shows  that  the  second 

nearest  neighbors  have  effects  of  the  order  of  4  *E  *  0.06  EQ  and  that  the 

_2 

third  nearest  neighbors  have  effects  of  the  order  of  6  EQ  *  0.03  EQ  for  a 
given  value  of  a/xQ  <  1.  This  means  that  the  inter-bubble  (and  also  inter- 
cloud)  interaction  is  short-ranged  and  that  the  third  nearest  neighbors  and 
beyond  have  no  significant  influence.  Thus,  only  a  small  number  of  bubbles 
are  necessary  to  model  an  N-bubble  system  with  N  »  1.  Note  that  the 
influence  of  the  bubbles  still  vanishes  at  infinity,  differing  from  systems 
satisfying  periodic  boundary  conditions. 

The  theoretical  treatment  described  is  exactly  applicable  to  any  N. 
However,  in  the  remainder  of  this  section,  we  include  up  to  a  total  of  5 
mutually  interacting  bubbles.  The  mathematical  manipulations  Involved  are 
analogous  to  those  of  the  preceding  section,  resulting  in  series 
expressions  similar  to  (19),  (20),  (24),  and  (25).  As  before,  only  dipole 
moments,  not  higher  multipoles,  are  Induced.  Since  no  new  insight  is  to  be 
gained  by  examining  the  actual  expressions,  we  give  below  only  the 
results.  In  figure  4,  we  show  the  electric  field  and  drift  configurations 
of  a  system  with  three  plasma  depletions.  The  three  depletions  are  again 
modelled  by  cylinders  of  radius  a,  located  at  x  ■  -  2xQ,  0,  and  2xQ.  Only 
one  quadrant  is  shown.  The  neighboring  bubbles  are  separated  by  a  distance 
2xQ  ■  2.5a,  and  are  100%  depleted  (K  ■  0)  as  before.  The  "external” 
electric  field  is  E^  given  by  (6),  and  the  density  profile  is  piecewise 
constant.  The  solid  lines  represent  the  polarization  electric 

A  2 

field  without  E  .  The  quantity  c  (E  -  E  )  x  B  /B  is  then  the 

— —  — o  — 

instantaneous  drift  velocity  relative  to  the  distant  undisturbed  ionosphere 
and  is  represented  by  the  dashed  lines.  Although  the  three-bubble  system 
is  different  from  a  two-bubble  system  in  that  the  former  has  a  central 
bubble  about  which  the  system  is  symmetric,  the  general  features  of  the 
field  and  drift  configurations  are  similar  as  can  be  seen  by  comparing 
figures  2  and  4.  That  is,  the  field  lines  and  drift  lines  inside  and 
around  the  end  bubbles  of  the  three-bubble  system  (figure  4)  are  similar  to 
those  of  the  two-bubble  system  (figure  2)  because  the  dominant  Influence 
arises  from  the  nearest  neighbors.  The  central  bubble  in  the  three-bubble 
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A  drawing  of  the  electric  field  lines  (solid)  and  the  x  drift 
lines  (dashed),  showing  one  quadrant  of  a  three-bubble  system. 
The  bubbles  are  placed  at  x  ■  +  2x0  and  x  ■  0.  The  points  A,  B, 
and  C  are  Inside  the  bubbles  at  x  -  0,  x  ■  a,  and  x  -  2xQ+a.  The 
separation  distance  Is  2xQ  -  2.5a  and  K  -  0  (100%  depletion). 


wmi. 


system  Is  affected  by  two  neighboring  babbles  with  comparable  effects.  The 
field  is  markedly  reduced  in  the  region  nearest  a  neighboring  bubble  (point 
B  in  figure  4). 

Figure  4  also  shows  the  distortion  in  the  electric  field  which  renders 
the  system  non-steady-state,  as  in  the  two-bubble  case.  The  drift  velocity 
(dashed  lines)  has  an  east-west  (horizontal)  component  that  may  be  a 
significant  fraction  of  the  vertical  velocity.  For  example, 

|E*/ (E*-Eq) !  *  0.48  at  9  *  n/8  from  the  point  B  just  Inside  the  boundary  of 
the  central  bubble  and  I E*/ (E*— Eq) I  »  0.55  at  0  ■  n/8  from  the  point  x  » 
xQ-a  and  y  *  0  just  inside  the  boundary  of  the  side  bubble.  Thus,  in  the 
neighborhood  of  these  points,  the  bubble  elements  should  have  significant 
horizontal  drifts. 

In  addition  to  the  distortion  of  bubble  contours  resulting  from  the 
non-uniform  electric  field,  the  figure  also  shows  that  the  electrLc  field 
is  substantially  reduced  from  that  of  an  isolated  single  bubble.  This  fact 
is  Illustrated  in  figure  5  which  gives  =  I /V ^ 1  for  the  points  A,  B  and 
C  corresponding  to  x  =*  0,  x  *  a  and  x  =  2xQ  +  a,  all  just  InsLde  the  bubble 
surfaces.  Here,  Vy  is  the  vertical  drift  velocity  of  the  three-bubble 
system  relative  to  the  distant  ionosphere.  Comparing  figure  5  with  figure 
3,  we  note  that  the  qualitative  behaviour  of  the  rise  velocity  is  similar 
in  both  systems,  exhibiting  significant  reduction  from  that  of  a  single 
bubble  system.  However,  in  the  three  bubble  configuration,  the  electric 
field  inside  the  central  bubble  is  substantially  weaker  than  the  two 
neighboring  bubbles.  Thus,  the  central  bubble  has  the  slowest  rise 
velocity.  The  line  D  in  figure  5  gives  the  relative  vertical  rise  velocity 
of  a  five-bubble  system  calculated  at  x  ■  a_.  Curve  D,  to  be  compared  with 
curve  B,  shows  that  the  Influence  of  the  additional  bubble  on  the  field 
inside  the  central  bubble  is  small.  The  reduction  in  rise  velocities  is 
increased  as  the  number  of  bubbles  is  Increased.  However,  the  influence  of 
the  bubbles  beyond  the  third  nearest  neighbor  is  small.  In  figure  6,  we 
show  the  relative  vertical  velocity  R  *  V^/V^  38  a  function  of  N,  the 
number  of  bubbles,  for  several  values  of  xQ.  For  a  given  separation 
distance,  the  rise  velocity  decreases  with  increasing  N  and  levels  off  for 
N  >  3.  Thus,  a  three-bubble  system  describes  well  the  basic  morphology  of 
an  N-bubble  system. 
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.g.  5  The  ratio  Rj  of  the  vertical  velocity  Vy  of  a  three-bubble  system 
to  the  single-bubble  rise  velocity  V^,  plotted  versus  xQ/a.  The 
curves  A,  B,  and  C  correspond  to  the  points  A,  B,  and  C  In  figure 
4.  The  curve  0  (dashed)  corresponds  to  Rj  ■  Vy/V^  at  x  »  a  In  a 
five-bubble  system,  to  be  compared  with  the  curve  B. 


Fig.  6  R  ■  Vy/V^  versus  N,  Che  number  of  bubbles,  evaluated  Inside  at  x 
•x0-a  for  the  two-bubble  case  and  at  x  ■  a  for  all  others.  K  ■ 
0.  The  separation  distances  (2xQ)  are  (a)  2x0  ■  2.5a,  (b)  2x0  - 
3.2a,  (c)  2x  *  4.0a,  and  (d)  2x  ■  10a. 


Figures  3  and  3  show  that  the  vertical  drift  velocity  of  a  bubble  is  a 
sensitive  function  of  the  separation  distance  (2x0)  except  for  relativly 
large  values  (2xQ/a  >6).  In  the  context  of  the  Raylelgh-Taylor 
instability,  this  implies  that  the  electric  field  configuration  and  the 
drift  velocities  may  depend  sensitively  on  wavelengths  (2xQ). 

As  a  general  remark,  we  point  out  that  the  preceding  results  derived 
for  plasma  density  depletions  for  which  K  <  1  are  also  applicable  to  plasma 
density  enhancements  (clouds)  for  which  K  >  1.  Calculations  for  the  cloud 
case  show  that  the  electric  field  configurations  are  qualitatively  similar 
to  that  of  the  bubble  case  (figures  2  and  4).  In  particular,  the  electric 
field  inside  the  clouds  experiences  the  greatest  reduction  In  the  regions 
facing  the  neighboring  clouds.  However,  the  boundary  condition 

K(E  )  ■  (E  )  implies  that  (E.  ).  Is  smaller  than  (E  ).  by  a 

in  i  out  1  v  in  L  v  out  L 

factor  of  K_1  for  clouds.  Thus,  for  a  given  separation  distance  2xQ,  the 
relative  distortion  in  the  electric  field  lines  Inside  a  cloud  is  less 
pronounced  than  in  a  bubble. 

Finally,  in  figure  7,  we  have  plotted  the  vertical  drift  velocity 
versus  K  5  n^/^  for  a  two-bubble  (K  <  l)  and  two-cloud  (K  >  1)  system. 
The  velocity  is  calculated  at  x  *  +  (xQ  -  a)  (point  A  in  figure  2)  and  Is 
upward  for  bubbles  and  downward  for  clouds.  The  velocity  is  normalized  to 
(equation  (27))  for  each  value  of  K.  Note  that  the  point  with  K  *  1 
does  not  exist  for  each  line.  For  K  *  1,  the  ionosphere  Is  not  disturbed 
and  there  is  no  bubble  or  cloud  drifting  vertically.  This  is  born  out  by 
the  fact  that  the  drift  velocity  vanishes  for  a  single-bubble  (cloud)  and 
any  multi-bubble  (cloud)  system.  Clearly,  the  velocities  vanish 
differently  for  different  separation  distances.  Mathematically,  the  ratio 
of  the  vertical  drift  velocity  Vy  -  c(Ex*  -  E0)/BQ  to  the  single-bubble 
(cloud)  drift  velocity  has  the  limit 

"  t1  “  rar  ®2f(x  +  vy)J 

as  K  +  1.  The  circled  points  in  figure  7  correspond  to  the  absolute  value 
of  the  quantity  which  has  no  physical  meaning. 
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Fig.  7  R2  “  Vy/V^  versus*  R  for  two-bubble  (K  <  l)  and  two-cloud  (K  >  l) 
systems.  The  separation  distances  (2xQ)  are  (a)  2xQ  -  2.5a,  (b) 
2xQ  ■  4a,  and  (c)  2xQ  ■  10a. 
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IV.  SUMMARY  AND  DISCUSSION 

In  Che  preceding  sections,  we  have  solved  the  current  conservation 

equation  (1)  with  quasi-neutrality  for  two-bubble  and  three-bubble  systems 
using  a  dielectric  analogy.  These  two  configurations  include  the  dominant 
near-neighbor  interaction  and  can  model  the  essential  morphology  of  the 
multi-bubble  systems  described.  This  is  demonstrated  by  actually 
calculating  the  field  including  up  to  five  bubbles.  The  solutions  are 
exact,  satisfying  the  specified  boundary  conditions  on  all  the  multiple 

disconnected  boundary  surfaces  and  at  infinity. 

Equations  (19),  (20),  (24)  and  (25)  give  the  solution  of  the  two- 

bubble  system  as  a  superposition  of  image  line  dipole  moments.  Similar 
expressions  are  obtained  for  three-  and  N-bubble  (cloud)  systems  (N  >  3). 
In  all  cases,  the  interaction  is  dominated  by  the  nearest  neighbors  and  is 
sensitive  to  the  separation  distance  (2xQ)  between  bubbles.  An  important 
result  is  that  the  electric  field  inside  the  bubbles  is  generally  highly 
non-uniform  so  that  the  multi-bubble  and  multi-cloud  systems  are  not 
steady-state  configurations  even  with  piecewise  constant  density  profiles 

(figures  2  and  4).  Moreover,  the  electric  field  inside  the  bubbles  is 
significantly  weaker  than  that  in  a  one-bubble  system  so  that  the  vertical 
drift  velocity  relative  to  the  undisturbed  ionosphere  is  slower  than  the 
one-bubble  case  (equation  (27)  and  figures  3  and  5).  For  moderately  small 
center-to-center  separation  distances,  2xQ  <  6a,  the  reduction  in  the  rise 
velocity  is  substantial.  This  implies  that  the  electric  field 

configuration  and  the  x  JJ  drift  velocity  may  depend  sensitively  on  the 
wavelengths  of  instabilities  causing  the  initial  density  fluctuations 
(e.g.,  the  Raylelgh-Taylor  instability).  In  particular,  the  bubbles  with 
smaller  wavelength-to-radius  ratio  would  rise  more  slowly  and  be  distorted 
more  strongly.  This  conclusion  may  be  particularly  applicable  to  the 
initial  linear  or  early-time  non-linear  stages  of  the  evolution.  it  has 
also  been  shown  that  in  some  regions  inside  bubbles,  the  horizontal  drift 
velocity  may  be  comparable  to  the  vertical  drift  velocity  with 

|E*/(e*-Eq)|  as  large  as  1/2  (Sections  II  and  ill). 

A  corollary  that  follows  from  the  non-uniform  electric  field  in  the 
bubbles  is  that  the  polarization  induced  J2  x  drift  velocity  is  sheared 
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and  that  a  component  of  the  electric  field  parallel  to  the  local  density 
gradient  Is  developed  In  the  leading  edges  of  the  babbles.  b'or  the  simple 
piecewise  constant  density  profiles  used  for  our  analysis,  the  density 
gradient  Is  not  well  defined.  In  a  more  realistic  density  profile, 
however,  some  previous  model  calculations  indicate  that  the  electric  field 
component  parallel  to  the  density  gradient  has  stabilizing  Influences  on 
the  Raylelgh-Taylor  instability  [Guzdar  et  al.,  1982]  and  the  £  x  IJ  drift 
instability  [Perkins  and  Poles,  1975;  Huba  et  al.,  1982).  lore 
specifically,  these  calculations  show  that  velocity  shear  preferentially 
stabilizes  the  short  wavelength  modes.  In  light  of  these  results,  we 
suggest  that  the  bifurcation  behavior  of  bubbles  and  clouds  may  be 
inhibited  by  the  presence  of  nearby  bubbles  (clouds). 

The  emphasis  of  our  analysis  has  primarily  been  on  plasma  density 
depletions  (bubbles).  However,  plasma  density  enhancements  (clouds  and 
striation  fingers)  can  also  be  treated  in  a  similar  fashion  [Scannapleco 
and  Ossakow,  1976;  Scannapleco  et  al.,  1976;  Ossakow  and  Chaturvedl, 
1978].  If  a  cloud  exists  in  the  equatorial  or  low  latitude  regie-'  ,  the 
preceding  results  are  all  applicable  with  the  replacement  of  K  >  1  (<  <  1 
for  bubbles)  where  K  -  n^/^.  In  particular,  as  thi  backside  an  initial 
cloud  begins  to  bifurcate,  the  small  xQ/a  and  roa  l  I  N  results  may  be 
applicable.  One  point  to  note  is  that  the  electric  field  inside  a  two-  and 
three-cloud  system  with  a  piecewise  constant  density  profile  is  similar  to 
that  shown  in  figure  2  and  4  with  weaker  field  strength  in  the  regions 
facing  the  neighboring  clouds.  However,  the  boundary  condition  K(E^n)^  » 
(Eout)j.  Implies  that  (Eln) L  is  smaller  than  (Eout)L  by  a  factor  of  K_1  for 
clouds.  As  a  result,  the  distortion  in  the  electric  field  lines  inside  a 
cloud  is  less  pronounced  than  in  a  bubble. 
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P.O.  BOX  808 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  TECH  INFO  DE? 
OICY  ATTN  DOC  CON  FOR  L-389  R.  OTT 
OICY  ATTN  DOC  CON  FOR  L-31  R.  HAGER 
OICY  ATTN  DOC  CON  FOR  L-46  F.  SEWAR 

LOS  ALAMOS  NATIONAL  LABORATORY 

P.O.  BOX  1663 

LOS  ALAMOS,  NM  87545 

OLCY  ATTN  DOC  CON  FOR  .J.  WOLCOTT 
OICY  ATTN  DOC  CON  FOR  R.F.  TASCHEK 
OICY  ATTN  DOC  CON  FOR  E.  JONES 

OICY  ATTN  DOC  CON  FOR  J.  MALIK 

OICY  ATTN  DOC  CON  FOR  R.  JEFFRIES 

OICY  ATTN  DOC  CON  FOR  J.  2 INN 

OICY  ATTN  DOC  CON  FOR  P.  KEATON 

OICY  ATTN  DOC  CON  FOR  D.  WESTSRVELT 

SANDIA  LABORATORIES 
P.O.  BOX  5800 
ALBUQUERQUE,  HM  87115 

OICY  ATTN  DOC  CON  FOR  W.  BROWN 

OICY  ATTN  DOC  CON  FOR  A.  THORNBROUG 

"OICY  ATTN  DOC  CON  FOR  T.  WRIGHT 

OICY  ATTN  DOC  CON  FOR  D.  DAHLGREN 

OICY  ATTN  DOC  CON  FOR  3141 
OICY  ATTN  DOC  CON  FOR  SPACE  PROJECT 

SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  B.  MURPHEY 
OICY  ATTN  DOC  CON  FOR  T.  COOK 

OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  D.C.  20545 

OICY  ATTN  DOC  CON  DR.  YO  SONG 

OTHER  GOVERNMENT 

DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON,  D.C.  20234 

(ALL  CORRESs  ATTN  SEC  OFFICER  FOR) 
OICY  ATTN  R.  MOORE 


/ 
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INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMUNICATIONS  &  INFO  ADMIN 
BOULDER,  CO  80303 

OICY  ATTN  A,  JEAN  (UNCLASS  ONLY) 

OICY  ATTN  W.  UTLAUT 
01CY  ATTN  D.  CROMBIE 
OICY  ATTN  L.  BERRY 

NATIONAL  OCEANIC  &  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 
OICY  ATTN  R.  GRUBB 
OICY  ATTN  AERONOMY  LAB  G.  REID 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
OICY  ATTN  I.  GARFUNKEL 
OICY  ATTN  T.  SALMI 
OICY  ATTN  V.  JOSEPHSON 
OICY  ATTN  S.  BOWER 
OICY  ATTN  N.  STOCKWELL 
OICY  ATTN  D.  OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

OICY  ATTN  RADIO  SCIENCES 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 

P.O.  BOX  983 
BERKELEY,  CA  94701 
OICY  ATTN  J.  WORKMAN 
OICY  ATTN  C.  PRETTIE 

BOEING  COMPANY,  THE 
P.O.  BOX  3707 
SEATTLE,  WA  98124 

OICY  ATTN  G.  KEISTER 
OICY  ATTN  D.  MURRAY 
OICY  ATTN  G.  HALL 
OICY  ATTN  J.  KENNEY 

BROWN  ENGINEERING  COMPANY,  INC. 

CUMMINGS  RESEARCH  PARK 
HUNTSVILLE,  AL  35807 

OICY  ATTN  ROMEO  A.  DELIBERIS 

CALIFORNIA  AT  SAN  DIEGO,  UNIV  OF 
P.O.  BOX  6049 
SAN  DIEGO,  CA  92106 


CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 
OICY  ATTN  D.B.  COX 
OICY  ATTN  J.P.  GILMORE 

COMSAT  LABORATORIES 
LINTHICUM  ROAD 
CLARKSBURG,  HD  20734 
OICY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 
OICY  ATTN  D.T.  FARLEY,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 
OICY  ATTN  H.  LOGSTON 
OICY  ATTN  SECURITY  (PAUL  PHILLIPS) 

ESL,  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 
OICY  ATTN  J.  ROBERTS 
OICY  ATTN  JAMES  MARSHALL 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 

OICY  ATTN  M.H.  BORTNER  SPACE  SCI  LAR 

CENERAL  ELECTRIC  COMPANY 
P.O.  BOX  1122 
SYRACUSF ,  NY  13201 
OICY  ATTN  F.  REIBERT 

GENERAL  ELECTRIC  TECH  SERVICES  CO.,  INC 
HMES 

COURT  STREET 
SYRYCUSE,  NY’  13201 
OICY  ATTN  G.  MILLMAN 

GENERAL  RESEARCH  CORPORATION 
SANTA  BARBARA  DIVISION 
P.O.  BOX  6770 
SANTA  BARBARA,  CA  93111 
OICY  ATTN  JOHN  ISE,  JR. 

OICY  ATTN  JOEL  GARBARINO 
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GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 
01CY  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
01CY  A 'TN  TECHNICAL  LIBRARY 
OICY  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

GTE  SYLVANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DZV 
77  A  STREET 
NEEDHAM,  MA  02194 

OICY  ATTN  MARSHALL  CROSS 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

OICY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 

OICY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
400  ARMY-NAVY  DRIVE 
ARLINGTON,  VA  22202 
OICY  ATTN  J.M.  AEIN 
OICY  ATTN  ERNEST  BAUER 
OICY  ATTN  HANS  WOLFARD 
OICY  ATTN  JOEL  BENGSTON 

1NTL  TEL  &  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

OICY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

11011  TORREYANA  ROAD 
P.O.  BOX  85154 
SAN  DIEGO,  CA  92138 

OICY  ATTN  J.L.  SPERLING 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAURAL,  MO  20810 

OICY  ATTN  DOCUMENT  LIBRARIAN 
OICY  ATTN  THOMAS  POTEMRA 
OICY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
OICY  ATTN  T.  MEAGHER 


KAMAN  TEMPO-CENTER  FOR  ADVANCEO  STUD: 
816  STATE  STREET  (P.O  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
OICY  ATTN  DASIAC 
OICY  ATTN  TIM  STEPHANS 
OICY  ATTN  WARREN  S.  KNAPP 
OICY  ATTN  WILLIAM  MCNAMARA 
OICY  ATTN  B.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 
OICY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
OICY  ATTN  DEPT  60-12 
OICY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  6  SPACE  CO.,  INC. 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

OICY  ATTN  MARTIN  WALT  DEPT  52-12 
OICY  ATTN  W.L.  IMHOF  DEPT  52-12 
OICY  ATTN  RICHARD  G.  JOHNSON  DEPT 
OICY  ATTN  J.B.  CLADIS  DEPT  52-12 

LOCKHEED  MISSILE  (.  SPACE  CO.,  INC. 
HUNTSVILLE  RESEARCH  6  ENGR.  CTR- 
4800  BRADFORD  DRIVE 
HUNTSVILLE,  AL  35807 
ATTN  DALF.  H.  DIVIS 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 

OICY  ATTN  R.  HEFFNER 

M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 

OICY  ATTN  DAVID  M.  TOWLE 
OICY  ATTN  P.  WALDRON 
OICY  ATTN  L.  LOUGH LIN 
OICY  ATTN  D.  CLARK 

MCDONNEL  DOUCLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
OICY  ATTN  N.  HARRIS 
OICY  ATrN  J.  MOULE 
OICY  ATTN  GEORGE  MROZ 
OICY  ATTN  W.  OLSON 
OICY  ATTN  R.W.  HALPRIN 
OICY  ATTN  TECHNICAL  LIBRARY  SKRVIC? 
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MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
OICY  ATTN  P.  FISCHER 
01CY  ATTN  W.F.  CREVIER 
01CY  ATTN  STEVEN  L.  GUTSCHE 
01CY  ATTN  D.  SAPPENFIELD 
01CY  ATTN  R.  BOGUSCH 
01CY  ATTN  R.  HENDRICK 
01CY  ATTN  RALPH  KILB 
01CY  ATTN  DAVE  SOWLE 
01CY  ATTN  F.  FAJEN 
01CY  ATTN  M.  SCHEIBE 
01CY  ATTN  CONRAD  L.  LONGMIRE 
01CY  ATTN  WARREN  A.  SCHLUETER 

MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

01CY  ATTN  JOHN  MORGANSTERN 
01CY  ATTN  G.  HARDING 
01CY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
01CY  ATTN  W.  HALL 
01CY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP. 

1456  CLOVERFIELD  BLVD. 

SANTA  MONICA,  CA  90404 
01CY  ATTN  E.C.  FIELD,  JR. 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

01CY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

442  MARRETT  ROAD 
LEXINGTON,  MA  02173 

OICY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 

OICY  ATTN  E.J.  FREMOUW 


PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

R  &  D  ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
OICY  ATTN  FORREST  GILMORE 
OICY  ATTN  BRYAN  GABBARD 
OICY  ATTN  WILLIAM  B.  WRIGHT,  JR. 
OICY  ATTN  ROBERT  F.  LELEVIER 
OICY  ATTN  WILLIAM  J.  KARZAS 
OICY  ATTN  H.  ORY 
OICY  ATTN  C.  MACDONALD 
OICY  ATTN  R.  TURCO 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
OICY  ATTN  CULLEN  CRAIN 
OICY  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

OICY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  END  AVENUE 
NEW  YORK,  NY  10023 

OICY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS,  INC. 

P.O.  BOX  2351 
LA  JOLLA,  CA  92033 

OICY  ATTN  LEWIS  M.  LINSON 
OICY  ATTN  DANIEL  A.  HAML1U 
OICY  ATTN  E.  FRIEMAN 
OICY  ATTN  E.A.  STRAKtK 
OICY  ATTN  CURTIS  A.  SMITH 
OICY  ATTN  JACK  MCDOUGALL 

SCIENCE  APPLICATIONS,  INC 
1710  COODRIDGE  DR. 

MCLEAN,  VA  22102 
ATTN:  J.  COCKAYNE 


SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 

OICY  ATTN  DONALD  NEILSON 
01CY  ATTN  ALAN  BURNS 
OICY  ATTN  G.  SMITH 
OICY  ATTN  L.L.  COBB 
01CY  ATTN  DAVID  A.  JOHNSON 
01CY  ATTN  WALTER  G.  CHESNUT 
OICY  ATTN  CHARLES  L.  RINO 
01CY  ATTN  WALTER  JAYE 
OICY  ATTN  M.  BARON 
OICY  ATTN  RAY  L.  LEADABRAND 
OICY  ATTN  G.  CARPENTER 
OICY  ATTN  G.  PRICE 
OICY  ATTN  J.  PETERSON 
OICY  ATTN  R.  HAKE,  JR. 

OICY  ATTN  V.  CONZALES 
OICY  ATTN  D.  MCDANIEL 

STEWART  RADIANCE  LABORATORY 
UTAH  STATE  UNIVERSITY 
I  DE  ANGELO  DRIVE 
BEDFORD,  MA  01730 
OICY  ATTN  J.  ULWICK 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

OICY  ATTN  W.P.  BOQUIST 

TRW  DEFENSE  A  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
RF.DONDO  BEACH,  CA  90278 
OICY  ATTN  R.  K.  PLEBUC1I 
OICY  ATTN  S.  ALTSCHULER 
OICY  ATTN  D.  DEE 

VISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
OICY  ATTN  W.  REIDY 
OICY  ATTN  J.  CARPENTER 
OICY  ATTN  C.  HUMPHREY 
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IONOSPHERIC  MODELING  DISTRIBUTION  LIST 
(UNCLASSIFIED  ONLY) 

COPY  TO  EACH  OF  THE  FOLLOWING  PEOPLE  (UNLESS  OTHERWISE 


PLEASE  DISTRIBUTE  ONE 
NOTED) 

NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  2037S 
DR.  P.  MANCE  -  CODE  4101 
DR.  E.  SZUSZCZEWICZ  -  CODE  4187 
DR.  J.  GOODMAN  -  CODE  4180 
DR.  P.  RODRIGUEZ  -  CODE  4187 

A.F.  GEOPHYSICS  LABORATORY 
L.G.  HANSCOM  FIELD 
BEDFORD,  MA  01730 
DR.  T.  ELKINS 
DR.  W.  SWIDER 
MRS.  R.  SAGALYN 
DR.  J.M.  FORBES 
DR.  T.J.  KENESHEA 
DR.  J.  AARONS 
DR.  H.  CARLSON 
DR.  J.  JASPERSE 

CORNELL  UNIVERSITY 
ITHACA,  NY  14850 
DR.  W.E.  SWARTZ 
DR.  R.  SUDAN 
DR.  D.  FARLEY 
DR.  M.  KELLEY 

HARVARD  UNIVERSITY 
HARVARD  SQUARE 
CAMBRIDGE,  MA  02138 
DR.  M.B.  McELROY 
DR.  R.  LIHDZEN 

INSTITUTE  FOR  DEFENSE  ANALYSIS 
400  ARMY/NAVY  DRIVE 
ARLINGTON,  VA  22202 
DR.  E.  BAUER 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
PLASMA  FUSION  CENTER 
LIBRARY,  NW16-262 
CAMBRIDGE,  MA  02139 


NASA 
Code  961 

GODDARD  SPACE  FLIGHT  CENTER 
CREENBELT,  MD  20771 
DR.  R.F.  BENSON 
DR.  K.  MAEDA 
Dr.  S.  Curtis 
Or.  M.  Dubin 

NATIONAL  TECHNICAL  INFORMATION  CENTER 
CAMERON  STATION 
ALEXANDRIA,  VA  22314 
12CY  ATTN  TC 

CO'CIANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20360 
DR.  T.  CZUBA 

COMMANDER 

NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIECO.  CA  92152 

MR.  R.  ROSE  -  CODE  5321 

NOAA 

DIRECTOR  OF  SPACE  AND  ENVIRONMENTAL 
LABORATORY 
BOULDER,  CO  80302 
DR.  A.  CLENN  JEAN 
DR.  G.W.  ADAMS 
DR.  D.N.  ANDERSON 
DR.  K.  DAVIES 
DR.  R.  F.  DONNELLY 

OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VA  22217 
DR.  G.  JOINER 

PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PA  16802 
DR.  J.S.  NISBET 
DR.  P.R.  ROHRBAUCH 
DR.  L.A.  CARPENTER 
DR.  M.  LEE 
DR.  R.  0I7ANY 
DR.  P.  BENNETT 
DR.  F.  KLF.VANS 
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PRINCETON  UNIVERSITY 
PLASMA  PHYSICS  LABORATORY 
PRINCETON,  NJ  08540 
DR.  F.  PERKINS 

SCIENCE  APPLICATIONS,  INC 
1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 
DR.  D.A.  HAMLIN 
DR.  L.  LINSON 
DR.  E.  FRIEMAN 

STANFORD  UNIVERSITY 
STANFORD,  CA  94305 
DR.  P.M.  BANKS 


UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 
DR.  N.  ZABUSKY 
DR.  M.  BIONDI 
DR.  E.  OVERMAN 

UTAH  STATE  UNIVERSITY 
4TH  AND  8TH  STREETS 
LOGAN,  UTAH  84322 
DR.  R.  HARRIS 
DR.  K.  BAKER 
DR.  R.  SCHUNK 


U.S.  ARMY  ABERDEEN  RESEARCH 
AND  DEVELOPMENT  CENTER 
BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN,  MD 

DR.  J.  HEIMERL 


UNIVERSITY  OF  CALIFORNIA, 
BERKELEY 

BERKELEY,  CA  94720 
DR.  M.  HUDSON 


UNIVERSITY  OF  CALIFORNIA 
LOS  ALAMOS  SCIENTIFIC  LABORATORY 
J-10,  M3-664 
LOS  ALAMOS,  NM  87545 
M.  PONGRATZ 
D.  SIMONS 
G.  BARAS CH 
L.  DUNCAN 
P.  BERNHARDT 

UNIVERSITY  OF  CALIFORNIA, 

LOS  ANGELES 
405  H1LLGARD  AVENUE 
LOS  ANGELES,  CA  90024 
DR.  F.V.  CORONITI 
DR.  C.  KENNEL 
DR.  A. Y.  WONG 


UNIVERSITY  OF  MARYLAND 
COLLEGE  PARK,  MD  20740 
DR.  K.  PAPADOPOULOS 
DR.  E.  OTT 
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